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■ Abstract. This review covers some of the developments in basic the- 

ory of astronomical masers over the past ten years. Topics included are 
d ' the effects of three dimensional geometry and polarization, with special 

emphasis on the differences between maser and non-maser radiation. 
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1. Maser Absorption Coefficient in Tliree Dimensions 



. The absorption coefficient k describes the coupling between particles and radia- 

I tion and serves as the foundation of maser theory. The appropriate expression, 

ly-^ . taking proper care of the Doppler matching between particles with thermal ve- 

I locity distribution and photons with a given wave vector k (aligned with the 

2;-axis) and frequency i', was developed by Litvak (1973). It can be written as 

• n\ ^ [ '^o(vk) , , 

7^. I/O 7 1 + /(vk)/Js 

O , 

f-j , where 
c/3 . 

Vk = {vx,Vy, = c[v - yQ\/vo), ^(vk) = ^ / ^(k')'^('^' - 1^0 + Vk • k')d^k'. 

^ \ Here /(k) is the intensity for wave vector k, is the saturation intensity and 

■ «^o(v) is the unsaturated absorption coefficient (proportional to the density of 

particles with thermal velocity v). Unfortunately, this expression was soon ig- 
nored and largely forgotten (except for a couple of papers by Bettwieser &: Kegel 
1974 and Bettwieser 1976) because shortly thereafter Goldreich &: Kwan (1974) 
introduced the much simpler expression 

(2) 



1 + JjJs ' 

where 

Koi/ = y" ^o(vk) dv^dvy , j 4^ ■ 

This expression gives the correct result for linear masers, where photon and 
particle motions are aligned, and because maser radiation is tightly beamed, 
Goldreich & Kwan reasoned that it should be adequate for all masers. It became 
the standard for all subsequent theory, including 3D geometries, even though it 
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was not derived from the proper equation [T|. It took many years until Neufeld 
(1992) recognized the internal inconsistency in this approach: in deriving the 
standard expression (|2|) the maser beaming angle is assumed to vanish, yet this 
very expression is used to solve the maser structure in any given geometry and 
derive the beaming angle for that geometry. 

Fortunately, the standard expression proved to be mostly adequate. Since 
maser radiation is tightly beamed, equation || can be handled with the aid of a 
series expansion in the beam width and the leading term in such a series repro- 
duces the standard expression 2 (Elitzur 1994). The standard theory provides 
the correct description of three dimensional masers and its results remain in- 
tact but only within the frequency core |x| ^ XsAud, where Avd is the Doppler 
width, X = {v—vq)//S.V£) and dimensionless parameter. For typical pump- 

ing schemes is ~ 2 in spherical masers, ~ 2.5-3 in disk masers and ~ 3-5 in 
cylindrical masers. For frequencies outside this core region, interaction with core 
rays that are slightly slanted to the direction of propagation suppresses photon 
production. Observed maser radiation is effectively confined to the core region, 
frequencies in the suppressed domain are essentially unobservable. In practice, 
suppression only affects extreme maser outbursts. Their profiles change in such 
a way that when fitted with a Gaussian they mimic line narrowing in proportion 
to (lni<o)~^/^, where Fq is the flux at line center, in contrast with the standard 
theory where such behavior is confined to unsaturated amplification. Such an 
inverse correlation between intensity and linewidth has been detected in a num- 
ber of II2O maser flares in star-forming regions (e.g. Mattila et al 1985, Rowland 
and Cohen 1986, Boboltz et al 1993, Liljestrom 1993). 

1.1. Linewidths in 3D Masers 

Maser radiation is tightly beamed, therefore the intensity 1^, angle-averaged 
intensity Jy and flux Fy (= / fil^dQ) are related via 

Fy = AnJy = ly^b (3) 

where ^Ib is the beaming solid angle. Saturated masers display two types of 
beaming (Elitzur, Hollenbach & McKee 1992). In matter-bounded masers, 
whose prototype is the filamentary maser, the beaming angle is controlled by 
the matter distribution and the maser observed size is equal to its physical 
size. The beaming angles of such masers are frequency independent, therefore 
the frequency profiles are the same for F^ and 1^. In amplification-bounded 
masers, whose prototype is the spherical maser, the beaming angle is controlled 
by the amplification process and the observed size is significantly smaller than 
the physical size. Because the amplification is strongest at line center the beam- 
ing is tightest there; the beaming angle increases with frequency shift from line 
center and the spectral shapes of F^ and are different from each other. 

With the standard expression for the absorption coefficient (|2|), it is easy 
to show that the flux of a saturated maser increases with length £ according 
to Fy oc KouJgi independent of the geometry (Elitzur 1990). Even with the 
full, proper equation |l|, the standard theory remains applicable at the core of 
the line so this result too is valid there. Therefore, in any saturated maser the 
frequency profile of the flux always obeys F^ oc kq^ oc exp(— x^), i.e., the flux 
spectral shape displays the full Doppler width Az^/j. On the other hand, the 
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intensity obeys oc Fi^/Qi^, therefore its spectral pprofilc will reflect also the 
frequency dependence of Qft. The result is I^, oc Kq^ oc exp(— ax^) where a is 
the dimensionality of the geometry; a is 1 for linear or filamentary masers, 2 for 
planar masers such as disks and 3 for fully 3D structures such as spheres. The 
width of the brightness spectral shape is thus 

A.= ^. (4) 

In particular, the intensity lincwidth of a saturated planar maser (the most likely 
geometry for shock induced masers) is 40% smaller than the Doppler width. 



2. Polarization 

Thermal radiation is generated in spontaneous decays, maser radiation in stim- 
ulated emission. This fundamental difference has profound implications, espe- 
cially for polarization. 

The stimulated emission process is the inverse of radiation absorption. Ab- 
sorption is a purely classical process, therefore the same applies also to stimu- 
lated emission. Since line radiation involves discrete energy states, the particle 
properties must be described with quantum mechanics. But the radiation wave- 
length is many orders of magnitude larger than particle dimensions so there is no 
need to quantize also the radiation field. The interaction of matter with maser 
radiation is adequately described with a hybrid, semi-classical approach (Litvak 
1970; Goldreich, Keeley k Kwan 1973, GKK hereafter): The radiation field is 
described by standard classical electromagnetic waves. The energy levels are 
eigenstates of the system Hamiltonian, treated by quantum theory. Interaction 
with the radiation field, treated as a perturbation, causes transitions between the 
energy levels. The transition rates for both absorption and stimulated emission 
are obtained from the product of the (quantum) matrix element of the transition 
dipole moment with the (classical) intensity of the radiation field. Since this is 
a complete description, not merely a classical analog, an important consequence 
is that there are no properties of the radiation generated in stimulated emission 
that are peculiar to the quantum theory; we must be able to fully deduce all of 
them from purely classical concepts as applied to propagating electromagnetic 
waves. 

In contrast, spontaneous emission is a purely quantum process. It has no 
classical analog since the initial state is devoid of radiation and thus cannot in- 
teract with any electromagnetic wave. Spontaneous decays do not occur even in 
standard treatments of quantum theory because the energy levels are stationary 
states of the system Hamiltonian, completely stable in the absence of external 
perturbations. This process occurs only when quantization of the electromag- 
netic field is taken into considerations, and can be interpreted as scattering off 
vacuum fluctuations. Spontaneous emission can be analyzed only in terms of 
the photon description of the radiation field. 

Induced Photons When stimulated emission is described in terms of photons, 
energy and momentum conservation imply that the induced photon has the 
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same frequency and direction as the parent photon. However, contrary to some 
widespread misconceptions, the induced photon does not have the same 

1. Phase: The argument of the oscihatory behavior of any wave is its phase 
(p = (j)Q + cut — li-r, where to is the angular frequency and k is the wave 
vector. An electromagnetic wave has a phase, a photon does not. The un- 
certainty principle leads to the relation AnAcp > 1 between phase and pho- 
ton number. The phase of a state with a well-defined number of photons 
is completely undetermined. Phases are meaningless when dealing with 
photon numbers, in particular they are irrelevant in spontaneous emission. 

2. Polarization: The induced photon polarization is not necessarily equal 
to that of the parent photon. Instead, it is determined by the change 
in magnetic quantum number m of the interacting particle. Am = 
transitions couple to photons linearly polarized along the quantization axis. 
Am = ±1 transitions couple to photons that are right- and left-circular ly 
polarized in the plane perpendicular to the quantization axis. Consider 
the interaction of a linearly polarized photon with particles in the upper 
level of a spin 1^0 transition. When the interacting particle is in the 
m = state it executes a Am = transition and the induced photon, too, 
is linearly polarized. But this is not the case when the particle is in one 
of the \m\ = 1 states. The linearly polarized photon, which can also be 
described as a coherent mixture of two circularly polarized photons, will 
now induce a |Am| = 1 transition and the induced photon is circularly 
polarized. Induced emission preserves polarization only when the magnetic 
transitions do not overlap. 

2.1. Polarization in Spontaneous Decays 

From Maxwell's equations, the electric field of an electromagnetic wave is always 
perpendicular to the propagation direction. This seemingly simple transverse 
condition is a rather peculiar constraint, hard to reconcile with the properties 
of the quantized particles that emit line radiation. Figure |l] shows the geometry 
for a Am = spontaneous decay. The quantization axis is denoted by B. The 
electric field generated in the transition, with an amplitude E^, is always parallel 
to this axis. The photons propagate in the direction marked by the double arrow 
at an angle 9 from B, the corresponding axis is denoted k. The component of 
the electric field along the axis parallel to the projection of B on the plane of 
the sky is E^^ = E^sinO. What about the component along the direction of 
propagation, E'fc? Is it E^ cos 9 as the geometry dictates? Or is it as required 
by the transverse condition? How can we reconcile these conflicting results? The 
answer is that we cannot as long as we apply classical reasoning. The resolution 
of this conflict is rooted in the quantum nature of spontaneous emission, which 
has no classical analog. Because of the uncertainty principle, only one component 
of any vector can be determined whenever the magnitude of that vector is known, 
the other two remain undetermined; recall the properties of angular momentum. 
The longitudinal component can be ignored in spontaneous emission — quantum 
mechanics can be counted on to take care of the transverse condition Ej. = Q 
and we can proceed directly to calculate the polarization. Denote by the 
intensity associated with the amplitude E^ [I^ oc Then the intensities 
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Figure 1. Polarization in Am = spontaneous emission. 

measured by a linear antenna oriented parallel and perpendicular to the B-axis 
are, respectively, Ty = sin^ 9 and Tj_ = 0. Therefore, the Stokes parameters of 
Am = spontaneous emission are / = Ty +Tj_ = I*^ sin^ 6 and Q = — T± = I, 
recovering the standard result of full linear polarization. 

2.2. Fully Resolved Zeeman Pattern; vb ^ 

When the magnetic field is sufficiently strong that the Zeeman shift vb exceeds 
the linewidth Ai^d, radiation is produced in pure Am transitions centered on 
the appropriate Zeeman frequencies. For the Am = transition we have just 
derived the polarization and it is straightforward to repeat these calculations for 
Am = ±1 spontaneous emission. The results are summarized in table |l] for the 
classical Zeeman pattern. In that case there are three spectral lines centered on 
1^0 + z^sAm (Am = 0, ±1), where uq is the line frequency in the absence of a 
magnetic field, with I^i^i') = + vb) = — vb)- Quantities listed in the 
first column are obtained for each transition from the product of the intensity 
heading the transition column with the appropriate trigonometric factor. The 
intensities that would be measured with right- and left-circular instrumental 
response are listed as T^- i and V = T^. — Ti. The parameter U vanishes for all 
transitions with this choice of axes. 

The results display the standard polarization properties of thermal radia- 
tion of fully resolved Zeeman components. And because each component can be 
considered an independent, isolated radiative transition that couples to a single 
sense of polarization, these results apply also to maser radiation even though 
they were derived for spontaneous emission. Indeed, these are the maser polar- 
ization properties derived by GKK, although from an entirely different approach. 
Thermal and maser polarizations are the same when the Zeeman pattern is fully 
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Table 1. Polarizations for fully resolved Zeeman pattern, 3> Ai'd 



resolved. The only difference between the two cases is the disparity between 
the vr and a maser intensities, reflecting their different growth rates (Elitzur 
1996). This disparity predicts a preponderance of cr-components. However, the 
evidence is mounting that individual vr-components are in fact never observed, 
a puzzle that currently has no explanation. 

2.3. Overlapping Zeeman Components; i'b ^ Ai//) 

The thermal and maser cases diverge now because the stimulated emission mixes 
the polarization components. Thermal emission is produced in spontaneous de- 
cays and must be considered in the photon picture; only intensities, i.e., photon 
numbers, are relevant. Maser radiation is generated in stimulated emission and 
its properties must be understood in terms of classical waves interacting with 
particles in quantized energy levels; both amplitudes and phases count. 



Thermal Radiation The various Am transitions produce spectral compo- 
nents with equal intensities but centered on frequencies slightly shifted from 
each other. The intensity of the Am = component is an even function of 
X, centered on x = 0. Introduce xb = vb/ Av^ <C 1, then I^{x ± xb) = I^{x). 
The three components are produced independent of each other and the overall 
radiation field is an incoherent superposition of them. Within each component 
the transverse condition is obeyed independently and the photons are polarized 
as described in table |l|. Adding up the Stokes parameters at every frequency 
across the line yields 



Q 

V 



1° + 1+ + r = 2/° 

g° + Q+ + Q- = [/° - V2(/+ + 1-) 



sm 



(Pi 

--^{XB sin 0)2 



+ + v- = (/+ - r 



I cos ( 



dl 

— xb cos I 

dx 



(5) 
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Figure 2. Geometry of the wave and quantization frames. 



The final expression in each case is the leading order result from a series ex- 
pansion in xb of I^{x) around /"^(x). The overall radiation field is polarized 
because it is the superposition of three polarized different intensities — because of 
the Zeeman shifts the intensities are slightly different at every given frequency. 
These differences are controlled by the parameter x b , and so is the polarization. 
When Xb — > 0, the polarization disappears. 



Overlapping Electromagnetic Waves At a given frequency and wave vec- 
tor, each Am (= 0,ibl) transition produces an electric vector with magnitude 



where the initial phase (pAm is random. The overall electric vector is E = 
J2 E^"^ and must obey the transverse condition E • k = 0. The geometrical 
setup is shown in figure |2[ Particle quantization is defined with respect to 
the x-y-z coordinate frame, with z the quantization axis. In this frame, each 
component of E is uniquely associated with a specific Am: the z-component 
couples only to Am = 0, so that = Ez, the x- and y-components couple 
only to Am = ±1, i.e., E^ = 2~^^'^{Ex ± iEy). The wave propagation is along 
the fc-axis, rotated by an angle 6 from the quantization axis in the x-z plane. 
Since the electric field is a proper vector, it can be decomposed in this frame 
too using straightforward, standard geometry. The transverse condition states 
that Ek = irrespective of the direction of propagation, and now we cannot 
rely on quantum considerations; this condition must be obeyed as a geometrical 
constraint on the three vector components of E (eq. ^). Given the amplitudes 
the condition Ek = becomes a relation among the phases 4'Am- For 
example, when [^''^l = \E~^\ = \E~\ the relation is 

(/)+ = -4>- , 4> = \(j)±\ = arccos (^2"^/^ cot 6^ , (7) 

where the meaningless overall phase is set through (pQ = (Elitzur 1991); one 
(T-component leads the vr-component by the phase difference and the other 
must trail by the exact same amount. Only waves launched with these phase 
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relations produce superpositions that are purely transverse so that they can be 
amplified by propagation in the inverted medium. The transverse components 
of these propagating waves are linearly polarized according to 

'? = 7 = -l + :rA^- (8) 
1 S sm u 

Whereas the thermal polarization arises from the superposition of different inten- 
sities, this one involves equal intensities but well defined phase relations among 
the amplitudes. The polarization arises because the independent constraints 
imposed by the particle interactions and the transverse condition must be rec- 
onciled simultaneously. Particle interactions (i.e., maser pumping) produce three 
independent fields corresponding to Am = 0, ±1. The transverse condition dic- 
tates that only two independent fields propagate in any given direction, the 
longitudinal combination of the original fields must vanish. The resulting phase 
relation, and polarization, refiect the correlation that must exist to eliminate 
the longitudinal component of E. The linear polarization in eq. ^ depends only 
on propagation angle, it is entirely independent xb- Indeed, the only assump- 
tion in its derivation was the existence of a quantization axis in the source — the 
physical process behind this axis was never specified, in principle it need not be 
a magnetic field. 

Equation |8| is immediately recognized as the polarization solution derived by 
GKK from an entirely different approach. Their assumption of equal pump rates 
for the different m-states is reflected in the equal l^^^™! taken here. The po- 
larization becomes unphysical {q > 1) for sin^ 9 < Y3. Only unpolarized maser 
radiation can propagate there since the interference dictated by the transverse 
condition cannot be obeyed for equal li?^™!. The polarization changes sign at 
sin^ 9 = ^/3, where it vanishes. At smaller (larger) angles q is positive (nega- 
tive), corresponding to polarization along (perpendicular to) the quantization 
axis. The transition between positive and negative q corresponds to a 90° flip 
in the polarization direction. Such flips are commonly observed in SiO masers; 
an example is shown in figure |3[ A natural explanation is a slight change in 
direction of the magnetic field, straddling the two sides of the transition angle 
9 = 55°. 
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Figure 4. The polarization vector 11 and the vector Kp that controls 
its radiative transfer in the space defined by the normalized Stokes 
parameters q = Q/I, u = U/I and v = V/I. 



It is important to note that the linear polarization will usually not exceed 
33% because only the limited range 35° < 9 < 45° gives \q\ > Y3. Propagation 
at 6* > 45° gives only l^l < Y3 — along the field when 45° < 9 < 55° and 
orthogonal to it when 9 > 55° . The direction of the magnetic axis is generally 
not known a-priori. It can be determined with certainty only when the linear 
polarization exceeds 33%, in which case the field projection on the plane of the 
sky must be parallel to the polarization. 

Maser Polarization The preceding discussion provides a geometric deriva- 
tion of the structure of polarization consistent with the fundamental physical 
processes that generate maser radiation. Equation |8| lists the only polarization 
consistent with the constraints that govern an interacting mixture of quantized 
particles and classical electromagnetic waves that have equal amplitudes in the 
quantization frame. In actuality, we cannot know these amplitudes beforehand 
and must derive them from a complete solution for the level populations coupled 
to the polarized radiation. The radiative transfer equation involves a matrix in 
the space of four Stokes parameter, presenting a rather complex problem. An el- 
egant geometrical interpretation was derived by Litvak (1975). The polarization 
structure of any electromagnetic wave is defined by the 3- vector of its normal- 
ized Stokes parameters 11 = (q, u, v) (figure ^). The off-diagonal elements of the 
radiative transfer matrix are ki oc ^/2{'n'i -|- nj~) — and Kc (x Ui — , where 
is the population of the magnetic m-state of the upper level. These elements 
can be combined to form another 3- vector Kp = {ki,0,Kc), then 

^ = [nx,.p]xn. (9) 

The effect of radiative transfer is to rotate the polarization vector of each individ- 
ual electromagnetic wave at the rotation velocity IIxKp; this velocity is different 
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for different waves and changes as 11 is rotating. Since waves are launched with 
arbitrary initial polarizations that subsequently rotate at different rates, the ra- 
diation field can be expected to remain unpolarized unless there is a stationary 
configuration whose polarization vector does not rotate. When such a config- 
uration exists, the polarization vectors are locked once they enter it and that 
becomes the polarization of the overall radiation field. 

Stationary polarization obviously occurs when 11 || Kp. However, since 
the level populations are affected by the interaction with the maser radiation, 
Kp itself varies and is affected by 11, therefore one must find a formalism to 
identify the stationary polarizations. It is straightforward to show that they 
are the eigenvectors of the radiative transfer matrix (Elitzur 1996). Two types 
of solutions exist for masers in a magnetic field. One type corresponds to the 
xb ^ I case, the other to ^ 1. The solution for the latter reproduces the 
GKK linear polarization accompanied by the circular polarization 

3cos^ 

This polarization arises when the populations of the magnetic sub-states of each 
level become equal to each other as a result of maser radiative interaction, re- 
quiring a unique cooperation between the particles and the radiation that they 
amplify. It is reached only in masers with Jg ^ 5, where S is the source func- 
tion, after the radiation has grown so that J/ Jg ^ xb- The condition q'^ + v'^ < 1 
constrains the propagation directions for polarized maser radiation. 

When the transition frequency varies, the Doppler width Aubi varies propor- 
tionately while the Zeeman splitting vb is unaffected. Therefore xb is inversely 
proportional to frequency, and the circular polarization decreases with the tran- 
sition frequency when all other properties remain fixed. Mcintosh, Predmore 
&: Patel (1994) find that SiO circular polarization indeed decreases when the 
rotation quantum number, and with it transition frequency, increases. The lin- 
ear polarization displays the opposite trend, increasing with rotation quantum 
number (Mcintosh & Predmore 1993). Since the solution linear polarization is 
independent of transition wavelength, this is the expected behavior in the pres- 
ence of Faraday depolarization, which is proportional to A^. The low rotation 
states are more severely affected because of their longer wavelengths and the 
linear polarization can be expected to decrease toward lower angular momenta, 
as observed. Although detailed calculations of Faraday rotation have yet to be 
performed for ^ 1, Mcintosh & Predmore find this to be the most plausible 
explanation of the data. 



2.4. Limitations and Outstanding Issues 

The theory presented here was developed for an idealized maser. The results 
depend in a crucial manner on the assumption of a constant direction for the 
quantization axis. They provide the maximal polarization that can be produced 
in a source that maintains a uniform direction for the magnetic field. Any curva- 
ture in the field lines along the propagation direction results in 6 variations that 
destroy the phase coherence between the vr- and ci-components, reducing the 
degree of polarization. In particular, Alfven waves introduce ripples in the field 
lines that destroy the polarization whenever the Alfven wavelength is shorter 
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than the amphfication length. Similarly, velocity blending will reduce circular 
polarization (Sarma, Troland & Romney, 2001) and linear polarization may also 
be reduced by Faraday rotation. As a result, the information that can be ex- 
tracted from polarization alone is limited because the same polarization can be 
produced in a number of different ways. For example, a certain linear polariza- 
tion can be attributed either to the maximal polarization at an appropriate angle 
or to a higher degree of polarization that was degraded either by curvature in 
the field lines or by Faraday depolarization. 

Another fundamental assumption is that the only degeneracy of the maser 
levels involves their magnetic sub-states. When any of the levels includes ad- 
ditional degeneracy, so that magnetic sub-states of different levels overlap, the 
tight constraints responsible for the stationary solutions no longer apply and the 
polarization can be expected to disappear]]. Indeed, the energy levels of both 
H2O and methanol involve hyperfine degeneracy and both masers are generally 
only weakly polarized. Exceptions do exist, though, and H2O masers sometime 
display high polarization, notably during outbursts such as in Orion (e.g., Abra- 
ham & Vilas Boas 1994). This may involve the excitation of a single hyperfine 
component, in which case the general solutions derived here are applicable. 



Incompleteness The theory presented here is incomplete. The radiation field 
is an ensemble of waves launched with random polarizations. Subsequent maser 
amplification through particle interactions is accompanied by rotation of each 
polarization vector, and we have identified the stationary modes that do not ro- 
tate. However, we have not shown how the radiation field actually evolves into 
this solution, and this is considerably more difficult. Indeed, it is always sim- 
pler to identify the stationary limit of a statistical distribution than to demon- 
strate how this limit is actually attained. Demonstrations of the approach to 
Maxwellian of a particle velocity distribution or to Planckian of a photon dis- 
tribution are considerably more difficult than the derivation of either functional 
form. The evolution of such ensembles requires numerical simulations of the 
type frequently performed in studies of plasma and laboratory lasers. This ap- 
proach can be avoided in the analysis of thermal radiation polarization, where 
phases are meaningless. But the essence of the maser polarization solution is 
specific phase relations among waves generated in different Am transitions, and 
those cannot be captured by standard radiative transfer techniques — the very 
derivation of the radiative transfer equation from Maxwell's equations is predi- 
cated on the assumption of random phases for different waves (cf Litvak 1970, 
GKK). A full simulation of the ensemble evolution of interacting particles and 
waves is the only way to study maser polarization growth. Such simulations have 
not yet been attempted for astronomical maser radiation. In addition to their 
inherent significance for demonstrating the approach to stationary polarization, 
these simulations are essential for a complete analysis of Faraday depolarization 
when I'B Az/£). 



^By example, consider the imaginary limit in which the hyperfine splitting of the OH molecule 
vanishes and the four ground-state lines are blended into one. 
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